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Abstract 

Left right symmetric models (LRSM) are extensions of the standard model by an enlarged gauge 
group SU(2)l x SU(2)r x U(1)b-l where automatic inclusion of right handed fermions as SU(2)r 
doublets guarantees a natural seesaw origin of neutrino masses. Apart from the extended gauge 
symmetry, LRSM also has an in-built global discrete symmetry, called D-parity which ensures 
equal gauge couplings for left and right sectors. Motivated by the fact that global symmetries 
are expected to be explicitly broken by theories of quantum gravity, here we study the effects of 
such gravity or Planck scale physics on neutrino masses and mixings by introducing explict D- 
parity breaking Planck scale suppressed higher dimensional operators. Although such Planck scale 
suppressd operators have dimension at least six in generic LRSM, dimension five operators can 
also arise in the presence of additional scalar fields which can be naturally accommodated within 
SO(W) grand unified theory (GUT) multiplets. We show that, such corrections can give rise to 
significant changes in the predictions for neutrino mixing parameters from the ones predicted by 
tree level seesaw formula if the left right symmetry breaking scale is lower than 10 14 GeV. 
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I. INTRODUCTION 



Left-Right Symmetric Models (LRSM) [1] have been one of the most well-motivated ex- 
tension of Standard Model (SM) studied in great details for last few decades. Apart from 
explaining the origin of parity breaking in weak interactions spontaneously, LRSM can also 
explain the origin of tiny neutrino masses [2] naturally via seesaw mechanism [3, 4] without 
reference to very high scale physics such as grand unified theories (GUT). Supersymmetric 
versions of such models have several other motivations like protecting the scalar sector from 
quadratic divergences, providing a natural dark matter candidate among others. However, 
as studied previously in [5, 6], generic Supersymmetric Left-Right models are tightly con- 
strained from consistent cosmology as well as successful gauge coupling unification point of 
view and in quite a few cases these models do not give rise to successful unification and 
consistent cosmology simultaneously. Recently, non-supersymmetric versions of LRSM were 
studied in the context of of gauge coupling unification and consistent cosmology [7]. It 
was shown that minimal versions of LRSM can not give rise to unification and consistent 
cosmology simultaneously, but suitable extensions of these models can give rise to both of 
these desired properties and at the same time allowing the possibility of a low scale gauge 
symmetry. 

Spontaneous breaking of exact discrete symmetries like parity (which we shall denote as 
D-parity hereafter) has got cosmological implications since they lead to frustrated phase 
transitions leaving behind a network of domain walls (DW). These domain walls, if not 
removed will be in conflict with the observed Universe [8, 9]. It was pointed out [10, 11] 
that Planck scale suppressed non-renormalizable operators can be a source of domain wall 
instability. The main theme of these works were to assume exact parity symmetry at tree 
level and introduce explicit parity breaking terms of higher order. As pointed out in [10], 
any generic theories of quantum gravity should not respect global symmetries: both discrete 
and continous. Without worrying about the details of such symmetry breaking mechanism, 
our purpose is to study the effects of such terms which arise only in the form of higher 
dimensional operators. The role of such higher dimensional operators in destabilizing domain 
walls was studied in [5, 7]. Here we intend to study the effects of such operators on the 
neutrino sector, namely the neutrino mixing parameters. We find that in generic LRSM, 
such operators which affect neutrino paramaters can have dimension at least six and do not 



2 



affect the neutrino masses and mixings significantly. However, in the presence of additional 
scalar fields, dimension five operators can arise and can affect the neutrino mass and mixings 
significantly. In particular, we incorporate the presence of gauge singlet scalar field which 
can naturally fit inside several SO(10) GUT representations. As discussed in [7], such singlet 
extension of minimal LRSM also leads to domain wall disappearance which is not possible 
in the minimal versions. Here we study the effect of such higher dimensional operators on 
neutrino mixing parameters and find that the corrections can be very significant if the left 
right symmetry breaking scale is below 10 14 GeV. 

This paper is organized as follows. In section II we discuss minimal LRSM with Higgs 
triplets and discuss how tiny neutrino mass arise in this model. In section III we discuss the 
possible higher dimensional and explicit parity breaking operators which can affect neutrino 
masses. Then in section IV, we present our numerical analysis on the effects of higher 
dimensional operators on neutrino mixing parameters and finally conclude in section V. 
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FIG. 1: Variations of Ara^ and Aro^ (NH), Am^(IH) as a function of f g for v R = 10 10 GeV 
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FIG. 2: Variations of Am^ and Am| 1 (NH), Am| 3 (IH) as a function of f g for v R = 10 12 GeV 



II. NEUTRINO MASS IN LRSM 



The fermion content of minimal LRSM is 

Ql= \ ^ j ~(3,2,1,^), Q R = ^ U * j ~(3*,l,2,i), 

^ J ~ (1,2,1,-1), £ R =l VR j 
eL y y e fi 

Similarly, the Higgs content of the minimal LRSM is 

* = | ^ | ~ (1,2,2,0) 

^12 012 

A L =\ L/V L ~ 1,3,1,2, A R =l * /V * -1,1,3,2 
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FIG. 3: Variations of Am^ and Am^(NH), Am^(IH) as a function of f g for vr = 10 GeV 

where the numbers in brackets correspond to the quantum numbers with respect to the 
gauge group SU(3) C x SU(2)l x SU(2)r x U(1)b-l- In the symmetry breaking pattern, 
the scalar A# acquires a vacuum expectation value (vev) to break the gauge symmetry of 
LRSM into that of the standard model and then to U(l) of electromagnetism by the vev of 
Higgs bidoublet $: 



SU{2) L xSU{2) R xU(l) B - L (A^ SU(2) l xU(1) y ($| U(l) em 
The relevant Yukawa couplings which leads to small non-zero neutrino mass is given by 

= VijiiL^ijR + y'i/iL$£jR + h.c. 

+ fa (C C ia 2 A R £ jR + (R^ Lj) + h.c. (1) 

where $ = r2$*r2. In the above Yukawa Lagrangian, the indices i,j = 1,2,3 correspond 
to the three families of fermions. The Majorana Yukawa couplings / is same for both left 
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and right handed neutrinos because of left-right symmetry. These couplings / give rise to 
the Majorana mass terms of both left handed and right handed neutrinos after the triplet 
Higgs fields ^ acquire non-zero vev. These mass terms appear in the seesaw formula as 
discussed below. The resulting seesaw formula in this minimal model can be written as 



Here mm is the Dirac neutrino mass matrix defined as ttilr = Uij{^)- It should be noted that 
the Yukawa couplings t/ij in the definition of Dirac neutrino mass matrix are not the same as 
the ones introduced in the Yukawa Lagrangian (1), but the ones obtained at the electroweak 
scale after renormalization group evolution (RGE) effects are taken into account from the 
scale of left right symmetry breaking down to the electroweak scale. Such RGE effects on 
neutrino parameters for type I and type II seesaw models have been studied in [12] and 
[13] respectively. However, in our present work we do not attempt to perform a systematic 
RGE study of neutrino parameters in LRSM. To simplify our analysis, we assume that the 
RGE effects on the neutrino Yukawa couplings from the left right symmetry scale to the 
electroweak scale do not diminish the effects of higher dimensional Planck scale suppressed 
operators on neutrino mass in LRSM (as we discuss in next section). 

In LRSM with Higgs triplets, Mrr can be expressed as Mrr = vrJr with vr being the 
vev of the right handed triplet Higgs field Ar imparting Majorana masses to the right- 
handed neutrinos and J'r is the corresponding Yukawa coupling. The first term Tnj^ L in 
equation (2) is due to the vev of SU{2) L Higgs triplet. Thus, = /lVl an d M RR = JrVr, 
where Vl^r denote the vev's and J'l,r are symmetric 3x3 matrices. The left- right symmetry 
demands /r = /l = /• The induced vev for the left-handed triplet vl can be shown for 
generic LRSM to be 

v L = 7 

vr 

with My/ ~ 80.4 GeV being the weak boson mass such that 

\v L \ « M w « \v R \ 

In general 7 is a function of various couplings in the scalar potential of generic LRSM and 
without any fine tuning 7 is expected to be of the order unity (7 ~ 1). The seesaw formula 



m LL = m£ 7 L + m LL 



(2) 



where the usual type I seesaw formula is given by the expression, 




-m L RM RR ml R 



(3) 
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in equation (2) can now be expressed as 

m LL = ^(M w /v R ) 2 M RR - m LR M^ml R (4) 



v R = 10 10 GeV 
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FIG. 4: Variations of sin 2 #i2 and sin 2 #23 as a function of f g for v R = 10 10 GeV 



III. HIGHER DIMENSIONAL OPERATORS IN LRSM 

In the minimal LRSM discussed above, the next to leading order terms contributing to 
neutrino masses can be written as 

C NR = f 9 d T iR C ia 2 k R l jR ^^ + R^L (5) 

m pl 

where Mpi ~ 10 19 GeV is the Planck scale. Here f g i ^ f g n and hence D-parity is explicitly 
broken with the introduction of the higher dimensional operators above. Now, using the tree 




FIG. 5: Variations of sin 612 and sin 623 as a function of f g for vr = 10 GeV 

level Yukawa terms (1) as well as the Higher dimensional operators (5), the right handed 
neutrino mass matrix can be written as 



M RR = fv R + fgR 



v R 



The first term on the right hand side of equation (2) takes the form 



™>LL = f V L + fgL 



V 



3 



M 2 Pl 



=► m£i = ^Mw/vRfifvR + f gL $j£) 

m pi 

=► mii = 1 (M w /vr)\M rr + f gL ^ - f gR ^) 

Thus there arises two additional terms in the seesaw formula after the higher dimensional 
operators are taken into account. These two terms are proportional to v\j '(vrM pi ) and 
vr/M pi respectively. We check that neither of these two correction terms can change the 
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FIG. 6: Variations of sin 612 and sin 623 as a function of f g for vr = 10 GeV 

predictions of neutrino parameters from the ones predicted by the tree level seesaw formula 
(4). This is obviously because of the l/M Pl supression in both the terms which is almost 
negligible compared to the tree level neutrino mass terms. 

Now, let us consider the presence of an additional gauge singlet field a in LRSM. Since 
a singlet like <r(l, 1, 1,0) can naturally fit inside several 5*0(10) representations, we assume 
the vev of this singlet field to be of order (a) ~ Mqut ~ 10 16 GeV. In the presence of such 
a field, the non-leading terms contributing to neutrino masses can be of dimension five as 
follows: 

C NR = f gL £j R C ia 2 A R £ jR -^- + R^L (6) 

Mpi 

Doing the same analysis as in the case of minimal LRSM, here Mrr is found to be 

(o)vr 



M RR = fVR + fgR- 
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The type II seesaw term m/ L becomes 

<l = l(M w /v R f(fv R + fj^p) 

Mpi 

=► mii = 1 (M w /v R ) 2 (M RR + (f gL - f gR ) { -^) 

Mpi 

Without losing any generality, we assume (f g L — f gR ) to be a Hermitian matrix of order one 
multiplied a numerical factor f g which decides the overall strength of the corrected term. In 
the next section, we study the variation of neutrino mixing parameters as a function of this 
numerical factor f g . 



IV. NUMERICAL ANALYSIS 



The latest global fit value for 3a range of neutrino oscillation parameters [14] are as 
follows: 

Amii = (7.00 - 8.09) x 10" 5 eV 2 

Am 2 31 (NH) = (2.27 - 2.69) x 10" 3 eV 2 

Am 2 3 (IH) = (2.24 - 2.65) x 10~ 3 eV 2 

sin 2 #i 2 = 0.27 - 0.34 

sin 2 # 23 = 0.34 - 0.67 

sin 2 # 13 = 0.016 - 0.030 (7) 

where NH and IH refers to normal and inverted hierarchy respectively. Unlike the tight 
constraints on the above parameters, the global fit 3a range for the value of Dirac CP phase 
5c p extends over entire — 2tt range. For illustrative purposes, here we take its value to be 
300 degrees (same as the central value given in [14]). 

For the purpose of our numerical analysis, we first fit the neutrino mass matrix ttill using 
the best fit global parameters mentioned above. For both normal and inverted hierarchical 
neutrino mass patterns, we consider extremal Majorana phases such that the mass eigenval- 
ues are either (m 1; m 2 , m 3 ) or (m 1; — m 2 , m 3 ) denoted by (+ + +) and (H h) respectively. 

We follow the same approach for numerical analysis as [15] where the variation of neutrino 
mixing parameters with respect to the dimensionless parameter 7 in the seesaw formula (4) 
was studied in details. 
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FIG. 7: Variations of sin 2 #i3 and sin 2 5cp as a function of f g for vr = 10 10 GeV 

After parametrizing the neutrino mass matrix for the tree level seesaw formula (4) using 
the global fit neutrino data, we introduce the correction term (6) to the seesaw formula. As 
discussed above, this correction term is of the form 

%T = l{M w /v R )\f gL - f gR )^ 

Mpi 

Here we assume (f g L — f g n) = f g O(l) where 0(1) is a Hermitian matrix of order one. 

We then vary the dimensionless parameter f g from 10 -5 to 1 and see the variations of 
neutrino mixing parameters. The results are shown in figures 1, 2, 3, 4, 5, 6, 7, 8 and 9 for 
three different values of left-right symmetry breaking scales vr = 10 10 , 10 12 , 10 14 GeV and 
both normal and inverted hierarchies as well as both types of extremal Majorana phases. 
As seen from the figures, the changes in the neutrino mixing parameters from the best fit 
values (corresponding to f g — in our case) become more and more significant as we go 
from vr = 10 14 GeV to vr = 10 10 GeV. In particular, for vr = 10 14 GeV, almost all the 
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FIG. 8: Variations of sin 2 #i3 and sin 2 5cp as a function of f g for vr = 10 12 GeV 

neutrino parameters lie within the 3a allowed range for all possible values of f g . Only Am^ 
goes outside the 3a range for f g > 0.1 for all the models and 6*13 deviates from the allowed 
range for f g > 0.35 for IH(+++) model. For vr = 10 12 , the mass squared differences lie 
within the allowed range only when f g < 0.01 whereas for vr = 10 10 GeV, they lie outside 
the 3a range for entire range of f g parameter under study. Similarly, the mixing angles are 
also found to lie within the allowed range only for some small range of parameter space for 
lower values of vr. 

V. RESULTS AND CONCLUSION 

We have studied the effects of higher dimensional Planck scale suppressed operators on 
neutrino masses and mixings in left right symmetric extension of standard models. These 
higher dimensional correction terms arise due to the fact that any theory of quantum gravity 
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FIG. 9: Variations of sin 2 #i3 and sin 2 5cp as a function of f g for vr = 10 14 GeV 

does not respect global symmetries: both continouous and discrete. Since left right symmet- 
ric models have an in-built discrete global symmetry called D-parity, it is generic to introduce 
explicit D-parity breaking terms suppressed by the scale of gravity or the Planck scale. We 
have shown that in the minimal LRSM, the order of such higher dimensional operators have 
dimension at least six and hence too small to affect neutrino masses and mixing. We then 
incorporate the presence of additional gauge singlet scalar field which allows dimension five 
Planck suppressed operators contributing to the neutrino mass matrix. Such gauge singlet 
field can be naturally fit within several SO (10) GUT multiplets. As discussed in our earlier 
work [7], these singlet scalar fields play a non-trivial role in destabilizing domain walls which 
arise in these models as a result of spontaneous D-parity breaking. 

Sticking to the issue of neutrino mass alone in the present work, we then fit the tree level 
neutrino mass matrix to the global best fit neutrino data. After doing this, we introduce the 
higher dimensional operators and see the variations in the neutrino mixing parameters with 
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the changes in the overall coupling strength of these opearators. We consider both normal 
as well as inverted hierarchies and two extremal Majorana phases in our work. Doing this 
exercise for three different left right symmetry breaking scales namely 10 14 , 10 12 , 10 10 GeV, 
we show that the effects of these operators can be very significant for those models with 
left right symmetry breaking scale below 10 14 GeV. It should be noted that the purpose 
of our study is not to rule out or disfavor any particular model, but to emphasize the fact 
that fitting the tree level seesaw formula with neutrino data is not enough in these models. 
The higher dimensional opearators which violate D-parity explicitly, can give rise to sizable 
contributions and hence have to be taken into account in generic left right symmetric models. 
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